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ABSTRACT 

In the synchrotron radiation model, the polarization property depends on both the 
configuration of the magnetic field and the geometry of the visible emitting region. 
Some peculiar behaviors in the X-ray afterglows of Swift gamma-ray bursts (GRBs), 
such as energetic flares and the plateau followed by a sharp drop, might by highly 
linearly-polarized because the outflows powering these behaviors may be Poynting- 
flux dominated. Furthermore, the broken-down of the symmetry of the visible emitting 
region may be hiding in current X-ray data and will give rise to interesting polarization 
signatures. In this work we focus on the polarization accompanying the very early 
sharp decline of GRB X-ray afterglows. We show that strong polarization evolution is 
possible in both the high latitude emission model and the dying central engine model 
which are used to interpret this sharp X-ray decline. It is thus not easy to efficiently 
probe the physical origin of the very early X-ray sharp decline with future polarimetry. 
Strong polarization evolution is also possible in the decline phase of X-ray fiares and 
in the shallow decline phase of X-ray light curves characterized by chromatic X-ray 
VS. Optical breaks. An XRT-like detector but with polarization capability on board 
a Swift-like satellite would be suitable to test our predictions. 

Key words: Gamma Rays: bursts — polarization — GRBs: jets and outflows — 
radiation mechanisms: nonthermal 



1 INTRODUCTION 

The successful launch of the Swift satellite in Nov 2004 opened a new window to reveal what has been happening in the early 
afterglow phase of gamma- ray bursts (GRBs). As summarized in Zhang et al. (2006) and Nousek et al. (2006), in a canonical 
Swift GRB X-ray afterglow lightcurve some interesting features are emerging (detected in a good fraction of but not all bursts), 
including a very early sharp decline (phase-I), a shallow decline or even plateau, (phase-II), preceding the conventional or 
so-called "normal" decay phase (phase-Ill & -IV), and the energetic X-ray flares (phase- V). Various modifications of the 
standard GRB afterglow model have been put forward to explain the observations (see Meszaros 2006; Piran & Fan 2007; 
Zhang 2007 for recent reviews). Although the underlying physical processes are not very clear, it is widely suggested that 
the X-ray Phases I, II and V may be relevant to prolonged activities of the central engine while the simultaneous optical 
afterglow is dominated by GRB forward shock emission (Fan & Wei 2005; Zhang et al. 2006; Ghisellini et al. 2007; Panaitescu 
2007). If so, the early time UV/optical and X-ray polarization behaviors should be independent and may be very different 
from each other. In some models, the outflows powering energetic flares and X-ray plateaus followed by a sharp drop are 
Poynting-flux dominated (e.g., Gao & Fan 2006; Giannios 2006; Troja et al. 2007), thus high-linear polarization of these 
synchrotron radiation photons is expected (Fan, Zhang & Proga 2005). 
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Except the configuration of the magnetic field, the geometry of the emitting area can influence the polarization property 
of the synchrotron radiation as well. The broken-down of the symmetry of the visible emitting region may be hiding in the 
peculiar X-ray data and can give rise to interesting polarization signatures. Motivated by this idea, in section 2 we calculate 
the polarization property of Phase-I and show that strong evolution is likely. In section 3, we argue that similar phenomena 
are also expected in the decline phase of the X-ray flares and possibly in the late time part of Phase-II, especially those 
associated with a chromatic break in optical band. We also discuss the detection prospect there. In section 4, we summarize 
our work with some discussions. 



2 THE POLARIZATION EVOLUTION ACCOMPANYING PHASE-I 

For phase-I, two leading interpretations are: (i) these X-ray photons are actually the high latitude emission of a main internal 
shock pulse (Fenimore et al. 1996; Kumar & Panaitescu 2000; Zhang et al. 2006) while the central engine has died; and (ii) 
these X-ray photons are powered by the weaker and weaker activity of the central engine, i.e., the sharp X-ray decline traces 
the activity of the dying central engine (Fan & Wei 2005; Piran & Fan 20070. The interpretation (i) has been widely accepted 
because some X-ray declines can be adjusted as oc — to)~^^'^'^^'' predicted by the high latitude emission model (Liang et al. 
2006; Yamazaki et al. 2006; Zhang, Liang & Zhang 2007a), where /3x is the X-ray spectral index, and to is a free parameter 
parameterizing the beginning of a certain internal shock pulsed. Such a test, however, is not conclusive. As shown in Sakamoto 
et al. (2007), for the very sharp X-ray decline, the data itself can not well constrain the function form {t ~ to)''' ■ In many 
events, k can be taken as a constant (for example, fc ~ 2) and thus is not relevant to f3x- Additionally, it is not easy to see why 
the GRB central engine turns off abruptly while the numerical simulation actually shows the contrast (MacFadyen et al. 2001; 
Zhang, Woosley & Heger 2007b). Interpretation (i) also fails to interpret some slowly decaying phase-I as in GRB 060614 
(|Butler fc Kocevski 20071 [Zhang et al. 2007a|IXu et al. 2008^ . As for interpretation (ii), there has been no strict observational 
test so far though it seems possible for any energetic burst phenomenon. 

We take into account both models in our calculation and examine the possibility to distinguish between these two models 
with future X-ray polarimetry. 



2.1 High latitude emission model 

If interpretation (i) is correct, i.e., the sharp X-ray decline is due to the high latitude emission of a main internal shock pulse, 
there would be very interesting polarization signals. As is known, GRB outflow is likely to be jetted. In this work we only 
consider the uniform jet model. Following Ghisellini & Lazzati (1999), we assume that the half-opening angle of the ejecta is 
6j, and that the line of sight (L.o.S.) makes an angle 9^ with respect to the jet's central axis (C.A., see Fig. la for illustration). 
The probability of observing the ejecta along its C.A. is vanishingly small since it corresponds a very small solid angle. For 
typical bright GRBs, the L.o.S. is likely to be within the cone, i.e., 9^ < 9y At early time, the high latitude emission is 
from zones around the L.o.S. satisfying !^ 9j — 0^, so the net polarization of the detected photons vanishes because of the 
symmetry; but at later time, the high latitude emission is from 9 > 9j — 9^ , the symmetry is broken and the observed net 
polarization is not zero any longer. 

In this model, the photons arrived at t were emitted from 

9 ^ {2ct/Rof^\ (1) 

where -Ro ~ -Rprompt is the radius of the main internal shock pulse and c is the speed of light. Below we adopt the simplified 
calculation of Ghisellini & Lazzati (1999). The notations used in Fig[l}3 are related to Rq, 9j and 9^ as R ^ JJosinSj and 
Til ~ Ro sin{9j — 9v). When p > R\, the ring is not a complete circle. The missing parts are in (0, V'l) a-nd {2-n — t/ji, 2-n) as 
shown in Fig[l] Approximately, we have (Ghisellini & Lazzati 1999) 

-iZ-i ^ i i " arcsin[^^(^j3^], for p > Rv, ^2) 
1 0, for p < Ri, 



^ These two interpretations can be also applied to the sharp decline of the X-ray flares. 

^ In this interpretation, the duration of the main prompt emission pulse should be in order of the timescale of whole burst otherwise 
the total high latitude emission of those early pulses decay too fast to account for the observation (see Fig. 2 of Fan & Wei 2005 for 
illustration). This is somewhat unusual in the standard internal shock model, in which the typical variability timescale is just in order 
of milliseconds. A possible solution of this puzzle is to assume that the prompt 7— ray emission are powered by some magnetic energy 
dissipation processes at a distance -Rprompt ~ 10^^ — 10^® cm to the central engine, so the emission duration is governed by the angular 
timescale ^ ^prompt /(2r2c) ~ 6(i?, prompt/10^"* cm)(r/50) ^ sec, where F is the bulk Lorentz factor of the outflow. 



Figure 1. (a) The high latitude emission model for the X-ray decline. The larger the angle with respect to the line of sight (L.o.S), the 
latter the emission arrives at us because of the geometry effect, (b) Sketch of the geometrical set-up used to compute the polarization 
signal (see also Ghisellini & Lazzati 1999). 



where p ~ Ro sin 9. 

In principle, the GRB outflow could be either baryon-rich or Poynting-flux dominated. For the former the magnetic fleld 
is generated in the shock front and is likely to be random. For the latter the magnetic field is coherent in a large scale. The 
resulting polarization hght curves are expected to be different in these two cases, as shown below. 

Random magnetic field. The net polarization of the emission from p — p + dp can be approximated by (See the 
Appendix for the derivation) 



^ n (2C'~l)sin^a [sin(27/>i)| 
^nPJ ~ ° 2<:2 sm^ a + {1 + cos^ a) 27r - ' ^ ' 

where Ho is the maximum polarization at one point, a is the viewing angle with respect to the fluid velocity in the outflow 
comoving frame (see Figure \KT\ for illustration), which is governed by cos a = {cos 9 — /3)/(l — l3 cos 9), and the parameter 
(^^ =< Bji > / < -Bx > denotes the level of anisotropy of the magnetic field distribution, where B± and By are magnetic fleld 
perpendicular and parallel to the normal of the shock plane. /3 = (1 — l/F^)^''^ is the velocity of the GRB ejecta in units of c. 

We note that in eq.^^ n(p) oc sin^ q/(1 + cos^ a) if ~ 0. For r{9j-9^) > 1, we have n(p) oc 2r^ {9j ~ 9^)^ /[I + r'^{9j - 
Sv)^] <S 1. The resulting net polarization can thus be neglected. But for 2(^^ sin^ a 2> 1 + cos^ a, the second term in the right 
side of eq.(|3} is — + 1. Significant net polarization is possible. 

As 9 > {9j — 9v), the observers just detected ~ (tt — '4^i)/tt part of the circle, the X-ray fiux will decline more steeply than 
^-(2+/3x) gjjjjj jg governed by 

i.^(i)ocI^^^t-(2+'3-). (4) 

TV 

We plot in Fig[2]the sharp X-ray decline and the corresponding polarization light curves. When the emission is from a 
circle ring, the net polarization is zero because of the symmetry. As 9 is getting larger and just part of the ring is visible, the 
observed polarization degree is not zero any longer until p^ — > 2RiR — Ri (i.e., tpi — > 7r/2). At the point of ipi = '''"/2, the net 
polarization becomes zero again according to equation (3), because the polarization direction abruptly changes by n/2. After 
this point, the polarization increases again. But the fiux may be too weak to have a significant polarization measurement. 

Homogenous magnetic field. In this case, the instantaneous stokes parameters (note that y = as the polarization 
is linear) and the polarization degree are given by (Granot & Konigl 2003): 

« } IS! ^ 

= no r2^-^i 7-7-777: ' (5) 



nZTV — pl r, l_y 



and 

n = ^^, (6) 

where 9p — (f> + arctan(i^ cotij)), y = V'^9^ and e = {1 ~\- l3x)/'2 (Granot & Konigl 2003). Again the polarization light curve 
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Figure 2. The sharp X-ray decline (upper panel) and the corresponding polarization light curves (lower panel) after a time-shift 
correction, in the case of the high latitude emission model. For this figure we assumed Ho = 60%, F = 50, 9j = O.f , 6v = 0.07, Rq = 10^^ 
cm and /3x = 1-15. 



has been plotted in FigE] which is similar to that of the random magnetic field case except an initial high polarization degree. 
This strongly suggests that the late time non-zero polarization is mainly contributed by the geometry effect, as in the case of 
a random magnetic field. If the symmetry has not been broken down, the late time net polarization for an ordered magnetized 
outflow disappears, as already found in Nakar, Piran & Waxman (2003). 



2.2 Dying central engine model 

Generally, in the dying central engine model there are two possibilities. (1) If the outflow powering the fast decaying X-ray 
tail emission always has a Lorentz factor F > l/{6j — dv), no strong polarization evolution is expected unless the physical 
composition of the outflow has changed a lot (e.g., from byron dominated to Foynting-flux dominated). We thus may be able 
to distinguish between the high latitude emission model and the dying central engine model with the future X-ray polarimetry. 
(2) The bulk Lorentz factor of the outflow drops so quickly as to satisfy F < l/{9j—9v) in a short period. When F < l/{9j—9v), 
the symmetry of the visible emitting area is broken and interesting polarization signature begins to exist. Since there is little 
information about the evolution of F with t, we here simply assume F « 150(t/10)"^ for 10 < i < 200 sec and 9j and 9v 
are constant. Again, we calculate the polarization evolution in the cases of random magnetic field and Homogenous magnetic 
field. But now the calculation should be due to the emission area within 6 ^ 1/F, rather than from a ring at 9 in the high 
latitude emission model. Therefore, we calculate the quantities of (Q, U, I) along the radius p and integrate them over p. Such 
calculations have been carried out by Ghisellini & Lazzati (1999) and Granot & Konigl (2003), respectively. Our numerical 
results are presented in Figure [S] One can see that polarization evolution is evident and shares some similarity with that in 
FiglH particularly in the case of random magnetic filed, though the polarization degree of high latitude emission model is 
larger than that of the dying central engine model. This is reasonable since in the former case the emission is from the ring 
while in the latter case the emission is from 9 ^ 1/F. In the case of homogenous magnetic filed, the evolution is small because 
of the significant polarization background (the degree is ~ 0.5) caused by the ordered magnetic field from the central source. 



3 PROSPECT OF DETECTING THE X-RAY POLARIZATION 

Our calculations in Section 2 are actually applicable to more phases in the GRB X-ray afterglow light curve. One is the decline 
phase of X-ray fiares (in Phase-V), which is also attributed to the high latitude emission model or the dying (re-)activity of the 
central engine. The other is the late phase of the X-ray flattening (in Phase-II). Recently Ghisellini et al. (2007) interpreted 
Phase-II, particularly those show chromatic breaks in X-ray and optical bands, as the prolonged emission of the central engine. 
The X-ray flattening ceases to exist when the bulk Lorentz factor of the outflow is so small that one sees the edge of the 
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Figure 3. The polarization light curves in the case of dying central engine model. The physical parameters are the same as that of FigO 
except we assume T = 150(t/10)-l for 10 < t < 200 sec. Please note that the dashed line is not the real polarization degree unless a 
factor 0.5 has been added. 



outflow. If these interpretations are correct, strong polarization evolution should be present at the end of Phase-II and in the 
decline of Phase- V. 

Measuring polarization is of growing interest in high energy astronomy. New technologies are being invented, and several 
polarimeter projects are under construction. In the X-ray band, the ongoing projects include XPE (Eisner et al. 1997), 
PLEXAS (MarsheU et al. 1998), POLAR (Produit et al. 2005), XPOL (Costa et al. 2006), POET (Polarimeters for Energetic 
Transients, see McConnel et al. 2007), and so on. For example, the planned POET includes two instruments: GRAPE (Gamma- 
Ray Polarimeter Experiment) and LEP (Low Energy Polarimeter, working in soft X-ray band). So it may be suitable to detect 
the polarization signatures predicted in the early X-ray afterglow of GRBs. An important issue is whether any of these detectors 
could perform a prompt slew to the GRBs localized by the 7— ray monitor. An ideal instrument would be an XRT-like detector 
with polarization capability on board a Swift-like GRB mission. 



4 SUMMARY & DISCUSSION 

Many surprises in GRB afterglows, mainly in the X-ray band, have been brought since the successful launch of Swift in 2004. 
Though not clearly understood yet, the (re-)activity of GRB central engine is likely to play an important role on producing 
afterglows, as speculated by Katz, Piran & Sari (1998). The different temporal behaviors in the UV/optical and X-ray 
afterglow favor that the peculiar X-ray emission traces the prolonged activity of the central engine while the optical afterglow 
is dominated by the forward shock emission UFan fc Wei 20051 IGhisellini et al. 2007|) . If true, the polarization behaviors of 
the early X-ray and optical afterglow may be independent and different. At such an early time, the linear polarization of 
the optical forward shock emission is expected to be very small (|Ghisellini fc Lazzati 1999]) unless the interstellar medium is 
magnetized ( [Granot fc Konigel 2003| ), because the Lorentz factor of the outflow is so large that the visible emitting region 
is well within the cone of the ejecta and is thus symmetric. But in the X-ray band, interesting polarization signatures may 
be present. For example. Fan et al. (2005) argued that some X-ray flares might have a high linear polarization because these 
new outflows are likely to be Poynting-flux dominated. Such an argument is also applicable to the X-ray plateau followed by 
a sharp drop as detected in GRB 070110. 

In this work, we focus on the polarization properties of the very early sharp decline of X-ray afterglows (Phase-I). 
In the high latitude emission model, the later the emission, the larger the angular 6 (see Figure [1] for details). So at late 
time, the symmetry of the emitting ring is broken when part of the emitting ring is out of the cone of the ejecta. Strong 
polarization evolution thus emerges (see Section 2.1). For comparison, results from the dying central engine model are a bit 
more complicated. If the outflow powering the fast decaying X-ray tail emission always has a Lorentz factor F > l/{9j — 9^), 
no strong polarization evolution is expected unless the physical composition of the outflow changed significantly. However, 
if the Lorentz factor drops rapidly to satisfy F < l/{9j — 9v), then polarization evolution is evident. In this scenario, if the 
magnetic field involved is shock-generated and is thus random, then strong polarization evolution is expected; while if the 
magnetic field involved is brought from the central engine and is thus large-scale ordered, then weak polarization evolution is 
expected due to the significant polarization background (the degree is ~ 0.5). One of our goals to calculate the polarization 
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of Phase-I is to see whether it's possible, via X-ray polarimetry, to distinguish between the high latitude emission model and 
the dying central engine model. This goal is not easy to fulfill in view of the above comparison. However, our calculation does 
support that strong polarization evolution is likely to accompany Phase-I, the decline phase of X-ray flares (Phase- V) and 
possibly the late part of Phase-II. Fruitful novel features are expected to appear in GRB X-ray afterglow polarimetry, though 
detection of them is beyond the capability of current monitors. We have witnessed the breakthrough made by Swift XRT and 
we are likely to experience it once more in the coming X-ray polarimetry era. One byproduct of the future X-ray polarimetry 
is to probe the cosmological birefringence effect that arises in some quantum gravity models (Fan, Wei & Xu 2007 and the 
references therein). 

Finally we'd like to point out that the optical flares in a few bursts (e.g., Boer et al. 2006) may also have a central engine 
origin (Wei, Yan & Fan 2006). Interesting polarization signatures are also expected. 
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Figure Al. The coordinate system to optimize the derivation (see also Sari 1999). 



APPENDIX A: DERIVATION OF EQUATION (3) IN THIS WORK 

Obviously, first we need to establish a suitable coordinate system to optimize the derivation. As stated in Sari (1999), at any 
point in the shock front there is a preferred direction, the radial direction, in which the fluid moves. We call this the parallel 
direction and choose the z-direction of the fluid local frame coordinate system to be in that direction. The two perpendicular 
directions {x,y) are assumed equivalent, i.e., the system is isotropic in the plane perpendicular to the direction of motion. We 
chose the a;-direction to be in the plane that contains the z-direction and the direction towards the observer n, as shown in 
Figure IXT] Suppose now that the magnetic field has spherical coordinates {6, ip) in that frame. A quite general description of 
the distribution of the magnetic field in such anisotropic system would be to allow different values of the magnetic field as 
function of the inclination from the preferred direction B = B{9) as well as a probability function for the magnetic field to 
be in each given inclination f{6). 

The relevant component of the magnetic field is that perpendicular to the observer i.e. Bsm{5), where S is the angle 
between the direction of the magnetic field and the observer. This will produce polarization Ho in the direction perpendicular 
both to the observer and to the magnetic field i.e. in the direction h x B. However, this polarization should be averaged due 
to contributions from magnetic fields oriented differently. On the other hand, because hght is a kind of transverse wave, its 
polarization should be perpendicular to its propagation. Therefore, we project the polarization to y (i.e., the component Q) 
and ri X y (i.e., the component U), respectively. 

Assuming that the emission is proportional to a certain power of the magnetic field, oc _B', the total polarization from a 
point- like region is theiH 




(Al) 



and 




(A2) 



and 




(A3) 



Since n — (sin a, 0, cos a) and B = (sin^cos(/), sin ^ sin 0, cos^), we have 
n = nx_B = (— cos a sin 9 sin 0, cos a sin 9 cos — sin a cos 9, sin a sin 9 sin 0) . 



With y = (0, 1, 0), we obtain 



cos rj = 



n-y 



cos a sin 9 cos — sin a cos 9 




[sin^ 9 sin^ (j) + (cos a sin 9 cos (f) — sin a cos 9)^]^^'^ 



where we have taken into account the fact that y, h X y and tl = h X B are all on the same plane which is perpendicular to n. 
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and 

sin 6 sin (b 

sin Tj ^ • 

[sin^ 6 sin^ + (cos a sin 9 cos — sin a cos 6)'^] ^/^ 

Prom the forms of cos rj and sin t] , we have 

^2 ^ 1 _ (cos a sin 9 cos — sin a cos 0) ^ — sin^ 9 sin^ 
sin^ 6 sin^ 4> + (cos a sin 9 cos — sin a cos 

and 

2 sin 6' sin (j!)(cos a sin 6^ cos — sin a cos 
sin^ 9 sin^ + (cos a sin ^ cos — sin a cos 

Now substitute equation (A5) into equation (A2), we get 



, ^ O ^ V V^WD Ct Dill (/ Olll CI V^WD 1/ I Dill 1/ Dill U/ /» j\ 

cos(277) = 2 cos^ r? - 1 = . Jo ' ( A4) 



sm(2r7) = 2 sm cos 77 = 2 ± \ ~ — ^ ■ ^ ■ (^5) 



= Ho / sin(2r?)[S(6i)sin5]V(6i)sin6ld0d6l 



n./ 
/ 



oc / sin 9 sin </)(cos a sin ^ cos — sin a cos 



/•TT /■27r 

/ sin 9 f{9)d9 / . 
Jo ^0 



sin 9 f{9)d9 / T{(p,a,9)d<f>. (A6) 


where !F(d),a,9) = , . ,^ . j^u^ have taken sin 5 = In x SI = ini. 

' [sm^ sin^ (i>+{cos a sin cos — sm CK cos (?)^] ^ ' ill 

To calculate the integration for J^{4>,a,9), we consider 

/■27r /•Ti' /'3-7r/2 /'27r 

/ T{(l),a,9)d(p=[j + + + ]T{4>,a,9)d(j}. (A7) 

In regions of ^ £ (0, 7r/2) and cji e (37r/2, 27r), sin(?i changes the sign but cos(/> does not. So in the above integral the 
contribution from the angle has been cancelled out by that from 2tt — (j). Similar results hold for the contribution of regions 
of € (7r/2, tt) and 4> € (tt, 37r/2). It is then straightforward to get 

p2n /'7r/2 p27v pTv c-Zi^ (2 

Uoc H<t>,a,6)d<f)=[(j + ) + (/ +/ )]J'i(l>,a,9)d(t> = 0. (A8) 

Jo Jo J:iTT/2 J-k/2 Jtt 

We thus have proved that because of the isotropy of the magnetic field in the {x,y) direction, the polarization of radiaiion 
emitted from a point-like region (after averaging on magnetic field orientation) must be in the direction perpendicular to the 
z-axis and to the observer, i.e., in the y direction. 

In the case of e = 2, both Q and / can be analytically integrated. With the relations J^" sin^ (pd(p — J^^ cos^ 0ci0 = %, 
we have 

PIT p2iv 

Q oc f (9) sin 9d9 / S^[(cosQ;sin^cos</) — sinacos^)^ — sin^ ^sin^ 

Jo Jo 



P7T 

TT sin^ a / - 
Jo 



f (9) sin 9[2 cos' 9- sin 9]d9, (A9) 



Note that < Sf >= B'^f{9) sin6>cos^ 9d9 and < Si >= B'^f{9) sin^ 9d9, we have 
< > 

Q oc 27rsin^Q(< By > - ^ ). 
Similarly, 

/"TT 

1=1 f{9) sin 9B^\U\'^d9 = 27r(sin^ a < Sjj > +(1 + cos^ a) < Bi > /2). 
Jo 



Therefore, the net polarization is given by (see also Gruzinov 1999 and Sari 1999) 

_Q _ sin^ a(< Bf > -< Si > /2) 

T ~ "sin^a < Bf^ > +(l + cos2a) < Bl > /2' 



(AlO) 



In the high latitude model for the X-ray sharp decline, the emission is from an angle a (or, equivalently p). As shown in 
Ghisellini & Lazzati (1999), it is convenient to write the polarization vector as a complex number P = Ppointe'^'"'' and integrate 
it between ipi and 2tt — ipi. Here is the position angle of the linear polarization in the observer frame. The polarization of 
a generic ring is thus 



9 



'''' -Ppoint/d6lp ^ ^ (2(^2 _ -^-j gjj^2 ^ I sin(2^i) 

r^""*' /de ^ ° sin^ a + (1 + cos^ a) 2-r - 2ipi 



n(p) = = no ,,,l^ , , i:r , (ah) 

'■/■I 

where =< By > / < ^1 >• 

Till here, we thus have proved equation (|3} in this work. 



